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In the paper we study the Breit-Wigner enhancement of dark matter (DM) annihilation con- 
sidering the kinetic decoupling in the evolution of DM freeze-out at the early universe. Since the 
DM temperature decreases much faster (as 1/R 2 ) after kinetic decoupling than that in kinetic equi- 
librium (as 1/R) we find the Breit-Wigner enhancement of DM annihilation rate after the kinetic 
decoupling will affect the DM relic density significantly. Focusing on the model parameters that try- 
ing to explain the anomalous cosmic positron/electron excesses observed by PAMEL A/Fermi/ ATIC 
we find the elastic scattering Xf — > Xf is not efficient to keep dark matter in kinetic equilibrium, 
and the kinetic decoupling temperature Tkd is comparable to the chemical decoupling temperature 
Tf ~ O(10)GeV . Fhe reduction of the relic density after Tkd is significant and leads to a limited 
enhancement factor ~ O(10 2 ). Therefore it is difficult to explain the anomalous positron/electron 
excesses in cosmic rays by DM annihilation and give the correct DM relic density simultaneously in 
the minimal Breit-Wigner enhancement model. 



I. INTRODUCTION 

The recent cosmic ray observations by 
PAMELA pj], ATIC0 and Fermi have all 
reported an excess of positrons and electrons 
from ~ 10 GeV up to ~ 1 TcV. These anomalies 
have stimulated a lot of interests, especially these 
excesses may be attributed to the signals of dark 
matter annihilation in the Galaxy. If these extra 
positrons/electrons are indeed from DM annihi- 
lation, it requires definite properties of DM. For 
example, DM should annihilate into lepton final 
states dominantly and should have a much larger 
annihilation cross section ((crv) ~ 10 23 cm 3 s _1 ) 
than the natural value ((crv) ~ 10 _26 cm 3 s _1 ) at 
freeze out The annihilation cross section at 

freeze out determines the DM relic density if DM 
is generated thermally at the early universe. 

In general, the DM annihilation cross section 
(av) depends on the averaged velocity of DM. For 
example, in the usual weakly interacting massive 
particle (WIMP) scenario, (av) can be expanded 
to a form of a+b(v 2 ) + 0(v A ) at the non-relativistic 
limit 0] • If the annihilation process is s-wave dom- 
inant, (av) is a constant. For the p-wave annihila- 
tion, (av) is proportional to (v 2 ). Therefore the 
DM annihilation by p-wave is suppressed today 
than the decoupling time since the WIMP usually 
has a velocity of v ~ 1CP 1 at the freeze-out epoch 
and cools when universe expands. The DM veloc- 
ity near the solar system is v ~ 10~ 3 , much smaller 
than that at the decoupling epoch. 

However, as indicated by the PAMELA, ATIC 
and Fermi data, we actually ask for a much larger 
annihilation cross section today to account for the 
excesses than that at the early universe. Contrary 



to the analysis before for the p-wave annihilation 
we require an annihilation form (av) depends on 
~ l/v n . This form leads to a large annihilation 
cross section today with low DM velocity and ex- 
plains the cosmic positron anomaly and relic den- 
sity simultaneously. Some mechanisms are soon 
proposed to achieve this aim after these results 
published, such as the Sommerfeld enhancement 
[1, Q and the Breit-Wigner enhancement [lfjl - flij . 

For the Sommerfeld enhancement, a new light 
mediator with mass of O(GeV) is introduced, and 
provide an enhancement factor of S ~ irax/v (ax 
is coupling constant between DM and mediator). 
For the Breit-Wigner enhancement, the DM an- 
nihilates via a narrow resonance, and an enhance 
factor of S ~ max[5, 7] _1 /O(10) can be obtained 
[3 (5, 7 are defined as S = (4m 2 - M 2 )/4m 2 and 
7 = Y/M respectively, where m is the mass of DM, 
M and Y are the mass and decay width of the res- 
onance respectively). One can achieve correct en- 
hancement factor S = (av) T=0 / (av) Tr ^ T (Tf is the 
temperature of chemical decoupling) by adjusting 
the parameters appropriately. 

It seems that the enhancement should not be 
important in the early universe when the veloc- 
ity of DM is ~ O(10 _1 ), and the enhancement 
factor is only S ~ 0(1). However, some recent 
studies showed that such effects are not negligible 
even at the freeze-out epoch [l5l - flll | , especially for 
the Sommerfeld enhancement. The Ref. [l7], EH 
pointed out that it may be difficult to achieve the 
required enhancement factor in the minimal Som- 
merfeld models considering the effect at the early 
universe. 

In this work, we will give a careful inspection 
on the Breit-Wigner mechanism at the DM freeze- 
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out process. For the Breit-Wigner mechanism, the 
DM annihilation continues after the chemical de- 
coupling until the DM velocity drops below the 
cut-off scale. Therefore the relic density is deter- 
mined by the cut-off scale related to 5 and 7 14 1. 



In the work we will show another important factor 
in determining the relic density, i.e. the kinetic 
decoupling process [21rJ23l ] . 

After the chemical decoupling at x ~ 20 (x rep- 
resents the temperature of the universe which is de- 
fined as x = m/T), the DM particle is still kept in 
kinetic equilibrium via the scattering with the hot 
bath. When such scattering is not efficient to keep 
DM in kinetic equilibrium, the DM momentum is 
red-shifted with the scale factor R, which leads to 
a rapid decrease of DM temperature as Tx ~ R~ 2 
rather than Tx ~ R~ x at the kinetic equilibrium 
epoch [2lU23j . Therefore, after the kinetic decou- 
pling (av) increases quickly and then reduces the 
abundances of DM more efficiently. Taking this 
effect into account we find the Breit-Wigner mech- 
anism is hard to provide a self-consistent explana- 
tion for both the DM relic density and the positron 
anomaly today. 

This paper is organized as following. In Section 
II, we briefly describe the Breit-Wigner enhance- 
ment mechanism at the DM freeze-out epoch. In 
Section III, we discuss the kinetic decoupling pro- 
cess. We will calculate the kinetic decoupling tem- 
perature and the DM relic density including such 
effect. In Section IV, we investigate the enhance- 
ment factor required by the cosmic positron mea- 
surements. We will study the parameter space and 
discuss whether there exists such parameters to ex- 
plain all the observations. Finally we give our con- 
clusions and discussions in Section V. 



II. THE BREIT-WIGNER 
ENHANCEMENT 

In Ref. [l2j, the DM annihilation process is as- 
sumed through XX —> R — >• //, w here R is a 
narrow resonance with mass M = ^jAm 2 (l - 5) 
and decay width V = M7 with |<5|,7 <C I. For a 
scalar resonance, the annihilation cross section is 
given as, 



I6tt 



1 



M 2 p i p i (<5 + v 2 /4) 2 + 7 



:B;Bf 



(1) 



where Pi and pi are defined as \/l — 4to 2 /M 2 and 
1 — 4m 2 /s respectively, s is given by s = (pi +P2) 2 , 
Bi and Bf denote the branching fractions of the 
resonance into initial and final states respectively, 
v is the relative velocity of two initial particles. For 



5 > 0, there exists an un-physical pole, but B-J 'pi 
is well defined. For simplicity, we parameterize the 
cross section as 11211 



(TV 



<5 2 +7 2 
(6 + z) 2 +' 



(2) 



Here a a = av\ T=0 means the cross section at zero 
temperature limit which is velocity independent, 
and is set as a free parameter in our work . z is 
defined in the form of s = 4m 2 (1 + z) which equals 
v 2 /4 in the non-relativistic limit. 

In order to calculate the DM relic density, it is 
necessary to solve the Boltzmann equation 
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\av)(Y 2 -Y 2 ) 



(3) 



where Y = hdm/s is the DM number density nor- 
malized by the entropy density s, A' is defined as 
A' = jj\ x=1 ■ The entropy density s(x) and the Uni- 
verse expansion rate H (x) of the universe are given 

by 
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where g*(gi) is the effective number of degrees of 
freedom for radiations (DM), and g*s is the ef- 
fective number of degrees of freedom defined by 
the entropy density. The (av) can be parameter- 
ized as (av) = o Q x~ n and the chemical decoupling 
temperature is obtained as 



xf ~ Ins — (n + l/2)ln(lne), 



(5) 



where e = c(c + 2)aA (c ~ 1 is a constant), A = 
X'a , and a = 0.145(<?i/g*) is defined in the form of 
e~ x at low temperature. The final Y 



ax 



3/2 



Y ei 

as x tends to 00 could be obtained approximately 
as Yoo ~ (n+ l)xy l+1 /A , and then the relic density 
n x h 2 = 2.74 x lO 8 ^^. 

In Ref. pjj, after parameterizing (av) for 5 > 0, 
the Boltzmann equation could be rewritten as, 



dY 
dx 



A 



S 2 + j 2 



x 2 (5 + &- 1 ) 2 +7 : 



■(r 2 -y e 2 ), (6) 



For the scalar resonance discussed above, <j 



<r is a combination of <5, 7 and other 
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parameters determined by the detailed model. It is in- 
deed a free parameter here. For more general discussions 
about the cross section formula of DM annihilation via 
s-channcl resonance, see Ref. Il9ll 
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where £ « l/\/2 is a constant (in fact, there is 



with 



an assumption here that x 



or DM stays in 



kinetic equilibrium until very low temperature in 
Eq. ([6])). For the Breit-Wigner enhancement, the 
freeze-out process begins at Xf ~ O(10) 2 , and con- 
tinues until the temperature of x& ~ max [5, 7] _1 
when the DM annihilation cross section does not 
increase with the universe cooling. The final value 
of Y is Yoo ~ Xf,/X. In the ordinary S-wave 
non-resonant annihilation scenario with (av) = 
constant, one could obtain Y^ ~ Xf /X , where 
x f0 ~ 20, A ~ A' x 10- 9 GeV^ 2 . Then the en- 
hancement factor is achieved as S ~ Xb/xfo — 
max[ < 5,7]- 1 /O(10) [3. The Breit-Wigner en- 
hancement has been used to explain the anoma- 
lous positron excesses which require an enhance- 
ment factor of ~ O(10 3 ). 



III. KINETIC DECOUPLING OF DM 
PARTICLES 

In the early universe, the DM production and 
annihilation processes XX ^ // are efficient to 
keep DM particles in chemical and kinetic equi- 
librium. After chemical decoupling at Tf DM 
may keep in kinetic equilibrium by momentum ex- 
change with the hot bath of the standard model 
particles via the t-channel scattering Xf Xf, 
until the temperature decreases to the kinetic de- 
coupling temperature Tkd- 

Before kinetic decoupling, the DM has the same 
temperature as the thermal bath. After kinetic de- 
coupling, the temperature of DM Tx decreases as 
1/R 2 , while the the temperature of thermal radi- 
ation still decrease as 1/R. So the Tx could be 
determined as I21H231 



T x — T, T x > Tkd 
T x = T 2 /T kd , T x < T h 



kd 



(7) 



Since Tx is different from T one can define a 
parameter x x related to DM temperature Tx as 



(8) 



where v is the most probable velocity of DM. The 
(av) is a function of x x which is given by (Tlj 



(av) 
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2 The x f could be achieved approximately by setting A 
A(<5 2 + 7 2 )/£ and n -v -2 in the Eq. ©. 
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where K\{x) and K 2 (x) are the modified Bessel 
functions of the first and second type respectively. 

After kinetic decoupling, the temperature of DM 
decreases rapidly, and the Breit-Wigner enhance- 
ment increases significantly. In the Fig. [TJ we 
show the enhancement factor of (av) / (av) \ x =20 f° r 
x kd = 50, 10 3 , 10 4 , 10 5 respectively. We also give 
the results in the limit of x kd = oo which denotes 
no kinetic decoupling. From Fig. [IJ we can see the 
(av) for x kd = 50 increases more quickly reaching 
the maximal value than the cases without kinetic 
decoupling. On the other hand, for a large value of 
x kd = 10 4 , 10 5 , such effects are not very obvious. 
These results could be understood easily from Eq. 
J2]) by assuming (av) ~ av\ z ^ v 2 roughly. When 
x > x kd and Vq ^> S,j, (av) increase as x 2 /x kd 



rather than x, and reaches ao more quickly for 
small x kd . 

The Fig. [5] shows the effects of kinetic decou- 
pling in the calculation of the relic density. Af- 
ter kinetic decoupling, the annihilation of DM be- 
comes more significantly, and reduce the relic den- 
sity more efficiently. If kinetic decoupling is very 
late x kd ^> Xf, for example x kd = 10 4 , such ef- 
fect is not very important compared with the case 
without kinetic decoupling. However, if the ki- 
netic decoupling occurs at nearly the same epoch 
as the chemical decoupling, the efficient annihila- 
tion would reduce DM relic density by about one 
order of magnitude. Therefore the kinetic decou- 
pling temperature T^d is a very important param- 
eter in the calculation of the DM relic density. 

The kinetic decoupling temperature Tkd can be 
determined using the method in Refs. [l8l |22|. If 
the momentum transfer rate drops below the ex- 
pansion rate, the DM decouples from the kinetic 
equilibrium with the radiation background. There- 
fore, the Tkd can be determined approximately by 
the relation of Tk(Tkd) = H(Tk d)- The momentum 
transfer rate is defined as Ref. UM, [231 



T 

T k ~ n r (av) s — 
m 



(12) 



where n r is the number density of masslcss 
fermions with n r = | • ^P-gfT 3 , (av)k is the 
thermally averaged cross section for the scattering 
process Xf — > Xf. Note that there exists a factor 
of T/m in the above formula, which reflects the 
approximate momentum transfer at each collision. 
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FIG. 1: The Breit-Wigner enhanced relative cross section (av) / (av) | x =20 as a function of x. The curves in the 
figures from left to right denote x kd — 50, 10 3 , 10 4 , 10 5 , oo respectively. The model parameters in the figures from 
left to right are set as 7 = S = 10~ 2 , 10 -5 , 10 -7 respectively. 



X kd =10 4 — 
X kd =1D 3 - 
X kd =50 - - 



FIG. 2: The evolution of DM abundance Vasa func- 
tion of x. The x kd are taken as 50, 10 3 , 10 4 , 00 respec- 
tively. 



Since the elastic scattering Xf — > Xf via t- 
channel is suppressed by the propagator of R with 
- A/ 2 ) 2 ~ 1/M 4 , the cross section of this pro- 
cess is much smaller than the annihilation cross 
section. The explicit formula of the cross section 
for Xf Xf depends on the model details. An 
approximate cross section of Xf —> Xf is related 
to XX // as 



a s a (5 2 +~/ 2 ) — , 



(13) 



T k (T kd ) = H(T kd ). Then we get 
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From above estimation, we can see the typical 
Tkd in the Breit-Wcigncr enhancement model is 
O(10)GeV, which is much larger than that in the 
ordinary WIMP model. For example, the T k d for 
neutralino in the SUSY model is only O(10)MeV 

mm. 

If T kd in Eq. flH]) is larger than the DM freeze- 
out temprature 3 Tf ~ rn/if ~ m/20, it means the 
clastic scattering becomes unimportant before the 
chemical decoupling. However, the DM particles 
are still kept in thermal equilibrium by the annihi- 
lation process // — > XX. Therefore T kd should be 
defined as min(Tf 1 T' kd ), where T' kd is determined 
by the elastic scattering as given in Eq. (|14p. 

For a more precise calculation, one need to de- 
rive the DM temperature T x (T) from the Boltz- 
mann equation 



L[f] = C[f] , 



(15) 



where the L and C are Liouvillc operator and col- 
lision operator for the scattering process respec- 
tively A general relation between the Tx and 
T has been provided by Ref. In our work, 

we still adopt the simple relation between the Tx 



where a s < 0(1) is a constant determined by the 

form of the interaction (for more details, see the 3 Here we define if as thc timc when „ dm = 
appendix). Then we can estimate Tkd by setting lOn^ (5/). 
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FIG. 3: Numerical illustration of the DM annihilation enhancement factor cr a / a nature without considering kinetic 
decoupling effect (left) and the corresponding x^d (right) on the 7 — 8 plane. 



and T as Eq. ([7]), and use the formulae in Ref. 
[HI to calculate the T' kd . We take a Z' model 
with m,DM = ITeV as an example, but our results 
can be extended to other models (for more details, 
see the appendix). From our calculations, we find 
that for the typical parameters used to explain the 
PAMELA/Fcrmi/ATIC results x kd = m/T kd is not 
far from if ~ O(10). To show this point explic- 
itly, we give the boost factor S = c /cr notere (left) 
and x kd (right) for different 5 and 7 in Fig. [3] Here 
"nature = 3x 10 -26 cm 3 s"' is the so called 'natural' 
value of DM annihilation cross section predicted by 
the WIMP models to generate correct relic density. 
In the left plot, we require each point in the pa- 
rameter space producing the correct relic density 
without kinetic decoupling effect, and determine 
the corresponding a . Then we use these <5, 7 and 
<t to calculate x kd . We find in the parameter space 
favored by the PAMELA/Fermi/ATIC results with 
S ~ O(10 3 ), x kd is similar as Xf. It means the ki- 
netic decoupling effect should be important in the 
early universe when determining the DM relic den- 
sity. Therefore it should be considered carefully in 
the explanation of the anomalous cosmic positron 
flux. 



IV. THE ENHANCEMENT FACTOR FOR 
ANOMALOUS POSITRON/ELECTRON 
FLUX 

In this section, we calculate the enhancement 
factor by the Brcit-Wigner resonance in the Galaxy 
today considering the kinetic decoupling. Here we 
define the enhancement factor as 



where a G denotes the {av) of DM with the most 
probable velocity v G ~ 10 -3 in the Galaxy. This 
definition is different from the earlier form s = 
G «l G na t U re E EH as the DM velocity is not zero 
today. We will see such difference is important. 

We give the numerical results of S in Fig. @] 
and Fig. [5] for the cases of 6 > and 5 < re- 
spectively. For each point in the two figures, <7q 
has been adjusted to produce the correct relic den- 
sity. The maximum value of S is only O(10 2 ) with 
5, 7 ~ O(10 -6 ) . From these results, we find the 
Breit-Wigner enhancement is difficult to provide 
large enough boost factor to explain the anoma- 
lous positron excesses after taking into account the 
kinetic decoupling effect. 

To check this result analytically, one would turn 
to the discussion in the last paragraph of Sec. |TT] 
fl2j . After the kinetic decoupling, the in Eq. 
([6]) should be modified by £xkd.x~ 2 - The DM an- 
nihilation would continue to the temperature of 
xi, ~ max[5, 7]" 5 • ^JXkd- One can also obtain an 
enhancement factor as S ~ Xb/xjo ~ max[S, 7]"' • 
y/xkd/xfo- It seems we could still achieve a re- 
quired boost factor by taking some smaller param- 
eters such as (5,7) ~ O(10 - ( 6 ~ 8 '). However, this 
is not the case. In fact, one can indeed obtain an 
arbitrary value of a / cr nature by setting the 6 and 7 
tiny enough as discussed above. But the factor of 
cr a /cr nature is different with that in the vanishing 
DM velocity limit. From the Eq. ^ we can see, 
for the parameters of6~7>O(10~ 5 ), these two 
factors are equal, because the (av) always reaches 
its maximum value a when the DM velocity de- 
creases to v„ ~ , 



S 



< max[5, 7] 2 . However, 
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is 



G 




s 

■100 

-80 



60 




FIG. 4: Numerical illustration of the enhance factor of 
S — <J G /o naturc on the 7 — 8 plane for the <5 > case. 




FIG. 5: The same as Fig. EJbut for 8 < 0. 



always smaller than er . 

To understand the maximum value of c G /<7 natuj , e 
and the corresponding parameters in Fig. 0] and 
Fig. [5J we can also use Eq. © as a roughly 
estimation. By setting z ~ v 2 G ~ 10 -6 and 

t7 /cr nature ~ max 7] _5 • y/xkd/xfo, we could ob- 
tain 5 

s _ max[(5,7]i ^x^ 
max[(5 + Vq, 7] 2 x/o 

From this rough estimation, we can see for the 
S > case, there actually exists a maximum value 
of S around max[S, 7] ~ w^. ~ O(10~ 6 ) as shown 
in the Fig. SJ On the other hand, for the 5 < 0, 
when 7 >• J + v 2 — > 0, the 5 might be larger 
than the case of S > 0. It means at the physical 



pole resonance, if the annihilations in the galaxy 
occur accurately with v 2 /A — 6—^0, the cross sec- 
tion could be very large. However, considering the 
dispersion of the DM velocity our numerical results 
show the enhancement factor can not be very large 
either. 



V. CONCLUSION AND DISCUSSION 

In this work, we study the Breit-Wigner en- 
hancement for DM annihilation taking the kinetic 
decoupling effect at the early universe into account. 
We find if the kinetic decoupling occurs at nearly 
the same epoch as the chemical decoupling, the 
DM annihilation process becomes very important 
and reduces the DM relic density significantly. Re- 
quiring the model gives correct relic density we find 
there is no parameter space that can give an anni- 
hilation cross section today large enough to explain 
the anomalous cosmic positron/electron excesses 
at PAMELA / ATIC /Fermi. 

The main point here is the elastic scatter- 
ing between DM and massless fermions Xf —> 
Xf is not efficient to maintain DM in ther- 
mal equilibrium. The kinetic decoupling occurs 
at high temperature ~ Tf. The DM temper- 
ature would decrease as ~ T 2 /Tf after kinetic 
decoupling, and reaches a very small value be- 
fore the structure formation. For typical WIMP 
such as neutralino, the typical damping mass is 
- 10- 6 M Q (m/100GeV)-^(T kd /30MeV)-% [H- 
Therefor in the Breit-Wigner enhancement with 
high Tkd , the damping mass might be much smaller 
than the usual cold DM model. This kind of 
DM model may predict tiny DM subhalo with 
M snJ < 1CT 6 M in the Galaxy. The realistic 
impact for the structure formation in the Breit- 
Wigner mechanism may need a careful study. This 
feature is possible to change the predictions for DM 
indirect detection. 

Finally we would like to point out that it is still 
possible to explain the anomalous cosmic positron 
excesses in some non-minimal Breit-Wigner mod- 
els. The ideal here is adding some new interaction 
process to keep DM in kinetic equilibrium till to a 
low temperature. For example, the DM is slepton 
t in the hidden sector [20]. It might interact with 
the hidden photon with a large coupling constant, 
or scatter with the standard model leptons by ex- 
changing hidden neutralino in resonance. The hid- 
den slepton annihilation to leptons could be en- 
hanced by a Z' resonance in the ?7(l)i j _i. model 
With this setting to enhance the scattering 
process, it is possible to obtain a low T^, and re- 
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cover all the discussions in the earlier works about 
the Breit-Wigner mechanism. 
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Appendix A: relation between cross sections of 
annihilation and elastic scattering processes 

In this appendix, we would give detailed discus- 
sions about the relation between the cross sections 
of annihilation XX — > ff and elastic scattering 
Xf -> Xf. 

We assume the effective interaction Lagrangian 
between two DM particles (X) and two leptons (f) 

as 



g A g B R XT x XfT f f 



(Al) 



where g A and gs are interaction couplings of XXR 
and ffR respectively, Tx and Tf are combines 
of Lorentz metrics determined by model, Rq is 
the propagator of resonance R. The cross section 
of annihilation process is (we neglect SM fermion 
mass m, here) 



1 



rg A g 2 B \R Qa \ 2 \M>\ 2 , (A2) 



da 

dfl 64tt 2 s Y s — 4m 2 

where we define the squared transition matrix el- 
ement is g A g B \Roa\ 2 \Ma\ 2 , s,t,u are usual Man- 
dclstam variables. For annihilation process in 
the non-relativistic limit, s is approximated as 
Am 2 + m 2 v 2 . So we can achieve a„ as 



1 



2 2 

9a9b 



r \M'\ 2 
■/LjUn, ( A3) 



327rm 2 M 4 (5 2 +7 2 ) 
For the cross section of elastic scattering, we find 
m . _ 1 9a9% f\M!> 2 



16ttto 2 M 4 J 4tt ~ dQ ' 
and the av s could be expressed by a as 

2 J\M>\ 2 dn 



(A5) 



In general, the \M' a \ 2 and \M' S \ 2 are expressed by 
the four vector momentums of four particles. We 



define p±, p 2 for two DM particles, and fa, fa for 
two fermions. We need only calculate either one 
of \M' a \ 2 and |Mg| 2 , and make some modifications 
to obtain the other one. In the calculation, we can 
neglect all the sub-leading terms which are pro- 
portional to m 2 , v 2 , and assumed the energy of 
fermion in the scattering process is u — 3T/2. 

For example, we can calculate the relation be- 
tween av s and a in a Z' model with Tx = = 
7**. The \M^\ 2 is given by 

\M' a \ 2 = i-32-[(fc 1 -fc 2 )m 2 + (p 1 -fc 1 )(p 2 -fc 2 ) 



+ (pi ■ k 2 )(p 2 ■ fa)] = 32m 4 



(A6) 



Then we can achieve \M'\ 2 = 8m 2 uj 2 and av 



h,(s 2 +i 2 )B 



Appendix B: calculation for the kinetic 
decoupling temperature 

In this appendix we show the calculation for the 
T kd in a Z 1 model. The detailed method is de- 
scribed in Ref. [23|, and can be extended to other 
models easily. 

In general, the DM temperature T x (T) can be 
derived by solving Boltzmann equation 



T v =T 



1 



Z n+2 



-exp 



[z]T[-(n + 2)- 1 ,z] 



(Bl) 

where z = (^;) n+2 - hi the low (high) temper- 
ature limit r^O (T — > oo), the T x has the same 
form T x -> T 2 /m (T x -> T) as Eq. Q. Then the 
kinetic decoupling temperature can be obtained as 



x / T-yO 



1% + 1 



(B2) 

The parameters a and n are defined as follows. 
One need to expand the amplitude at t = and 



2mui 



\M\ 2 = c n 



o | | 

TO 



(A4) The constant a is given by 



10 



(2tt)V 



1/2 



g f c n N n+3 — 



(B3) 



(B4) 



The N^ +3 for fermion (plus sign) and scalar (minus 
sign) are given by 



JV :: 



(l-p ± 2-")(n+l)!C(n + l) 



(B5) 



8 



where p + = 1 and p =0. For a Z' model with 
resonance mass mz' = M ~ 2m described as 



9a9b 
M 2 



(B6) 



we can achieve the amplitude at zero momentum 
transfer as g 2 A g 2 B \M' s \ 2 =0 /M i = ^fi-(^) 2 - Substi- 
tuting n = 2 and C2 = 9a9b /2 in the Eq. (|B2|) , we 
obtain T, . 



T,. 



1.326 



TO 



514 



S/C2M 



(B7) 



Here we assume the has the same couplings with 
the different leptons and sum the gj together. The 
thermal average annihilation cross section at low 
temperature can be written as <7o = C2/87tto 2 (<5 2 + 
7 2 ), the Eq. (|B8|) can be re- written as 



T, 



0.6 



g*vn 



g f a {8 2 +1 2 )M p: 



which is smaller than the result from Eq. 
a factor of 0(1). 



(B8) 
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